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Abstract 

In this paper, we introduce two discrete curvature flows, which are called u-fLows 
on two and three dimensional triangulated manifolds. For triangulated surface M, we 
introduce a new normalization of combinatorial Ricci flow (first introduced by Bennett 
Chow and Feng Luo [3]), aiming at evolving a order discrete Gauss curvature to a con¬ 
stant. When ax{M) < 0, we prove that the convergence of the flow is equivalent to the 
existence of constant u-curvature metric. We further get a necessary and sufficient 
combinatorial-topological-metric condition, which is a generalization of Thurston’s 
combinatorial-topological condition, for the existence of constant a-curvature metric. 
For triangulated 3-manifolds, we generalize the combinatorial Yamabe functional and 
combinatorial Yamabe problem introduced by the authors in to a-order. We also 
study the a-order flow carefully, aiming at evolving a order combinatorial scalar cur¬ 
vature, which is a generalization of Cooper and Rivin’s combinatorial scalar curvature, 
to a constant. 


1 Introduction 

Calculations in local coordinate charts are indispensable for the study of smooth mani¬ 
folds. However, an entirely different way came from Regge m- The basic procedure is 
to triangulate a manifold to simplexes, and finally obtain the geometrical and topological 
information of the manifold by studying the geometry in a simplicial complex and the 
combinatorial structure in the triangulation. On triangulated manifolds, the most general 
discrete metric seems to be piecewise linear metric which is defined on all edges. Discrete 
curvatures are determined by discrete metrics. Roughly speaking, discrete metrics are 
edge lengthes while discrete curvatures are angles between sub-simplexes and their com¬ 
binatorics. Besides the general piecewise linear metrics, there are discrete metrics dehned 
on all vertices. This type of metrics may be considered as a discrete conformal class. In his 
work on constructing hyperbolic metrics on 3-manifolds, Thurston m introduced circle 


1 


packing metric on a triangulated surface with prescribed intersection angles. For triangu¬ 
lated 3-manifolds, Cooper and Rivin introduced a sphere packing metric. We shall study 
circle (sphere) packing metrics and the corresponding discrete curvatures in the following. 

Suppose M is a closed surface with triangulation T = {V,E,F}, where V,E,F repre¬ 
sent the sets of vertices, edges and faces respectively, ^ [0, is a function evaluating 

each edge i ~ j a weight ^ij. The triple (M, T, <h) will be referred as a weighted trian¬ 
gulation of M in the following. All the vertices are supposed to be ordered one by one, 
marked by 1, • • • ,N, where N = |1/| is the number of vertices. Throughout this paper, all 
functions / : R —)• R will be regarded as column vectors in and /* is the value of / at 
i. And we use CiV) to denote the sets of functions defined in this way. 

Each map r : R —>■ (0, -|-oo) is called a circle packing metric. Given (M, T, ‘h), we can 
attach each edge {ij} a length h = + r'j + 2rirj cos ^ij. It is proved [l9] that the 

lengths {lij-,ljk-,hk} satisfies the triangle inequality for each face {ijk} G F, which ensures 
that the face {ijF} could be realized as a Euclidean triangle with lengths {hj^ljk^hk}- 
Suppose Oi^ is the inner angle of triangle {ijk} at the vertex i, the classical discrete Gauss 
curvature at the vertex i is defined as 

K, = 2tt- Y. di\ ( 1 . 1 ) 

{ijkjeF 

where the sum is taken over all the triangles with i as one of its vertices. Then 

N 

YKi = 27rx{M), 
i=l 

which may be considered as a discrete version of Gauss-Bonnet formula. Denote Kav = 
2'nx{M)/N. Notice that constant A'-curvature metric, i.e. a metric r with Ki = Kav for 
all i G V, does not always exist. In fact, for any proper subset I C R, let Fj be the 
subcomplex whose vertices are in I and let Lk{I) be the set of pairs (e, v) of an edge e and 
a vertex v such the following three conditioins: (1) the end points of e are not in /; (2) v 
is in I; (3) e and v form a triangle. Thurston |19j proved that the existence of constant 
A'-curvature metric is equivalent the following combinatorial-topological conditions 

2 vrx(M)M>_ ^ (7r-$(e))+27rx(A,), VI : 0 C / R R. (1.2) 

{e,v)£Lk{I) 

Moreover, the constant AT-curvature metric is unique, if exists, up to the scaling of r. 

Ghow and Luo [3] first established an intrinsic connection between Thurston’s circle 
packing and surface Ricci flow. They introduced a combinatorial Ricci flow 

F = -if... (1,3) 
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and its normalization 


QjV ' 

= {Kav - Ki)ri. (1.4) 

Then they proved that flow ()1.4p converges iff. the constant ii'-curvature metric exists, 
iff. Thurston’s combinatorial and topological conditions (II.2p are satisfled. 

Inspired by Chow and Luo’s work, the first author introduced a combinatorial Calabi 
flow 

^ (1.5) 

in [5] and proved similar results. The authors also studied the corresponding problem in 
hyperbolic background geometry using combinatorial Calabi flow [6]. 

The paper is organized as follows. In Section 2, we introduce the ath order Yamabe 
flow for surfaces and use it to study the corresponding constant and prescribing curvature 
problem. In Section 3, we study the constant a-curvature problem for triangulated 3- 
manifolds. And Section 4 is devoted to some useful lemmas used in the paper. 


2 a order combinatorial flow in two dimension 


As pointed by the authors in [9], there are two disadvantages of the classical definition of 
Ki. For one thing, classical discrete Gauss curvature does not perform so perfectly in that 
it is scaling invariant, i.e. if fj = Arj for some positive constant A, then Ki = Ki, which 
is different from the transformation of scalar curvature R\g = \~^Rg in smooth case. For 
another, classical discrete Gauss curvature can’t be used directly to approximate smooth 
Gauss curvature since it always tends to zero as the triangulation becomes finer and finer. 
To amend this flaw, we suggested a new definition of discrete Gauss curvature by dividing 
an “area element” Aj. 

We observed that the easiest form of “area element” may be Ai = vrr“, where a is any 
real number. For convenience, we omitted the coefficient vr in the following definition. 


Definition 2.1. Given a weighted triangulated surface (M, T, <h) with circle packing met¬ 
ric r : 1/ —>■ (0, +oo), the discrete Gauss curvature of order a (“a-curvature” for short) at 
the vertex i is defined to be 

( 2 . 1 ) 


Ra,i — 


rr^OL ’ 


where Ki is the classical discrete Gauss curvature defined as the angle deficit at i by dEH). 


For the classical discrete Gauss curvature K, [miiais] provide a complete description 
of the space of admissible LC-curvature {K = K{r)\r E see also Lemma 14.51 in 

this paper. Similarly, it’s interesting to know how to describe the space of admissible 
a-curvatures. 
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Remark 1. For the special case a = 2, Ai = vrr? is just the area of the disk packed at 
i. This case is especially interesting. It was shown [9] that rf is a suitable analogue of 
smooth Riemannian metric tensor g. Furthermore, the 2-curvature is very similar to the 
smooth Gauss curvature. On one hand, it has similar scaling law; On the other hand, it 
can be used to approximate the smooth Gauss curvature. 

One may ask if there are circle packing metrics whose corresponding a-curvatures 
are constants or have prescribing curvatures. Gurvature flow method gives an efficient 
approach to study these questions. The smooth 2-dimensional Yamabe flow, which is the 
same as Ricci flow, is dgjdt = (r — R)g. It is Ghow and Luo’s insight to consider r* as an 
analogue of smooth Riemannian metric tensor g, Ki as an analogue of smooth Riemannian 
curvature K = Rj2 and Kav as an analogue of smooth average curvature r. Thus the 
classical combinatorial Ricci flow (jl.4p comes into being. Inspired by their work and the 
idea of taking r? as a suitable analogue of smooth Riemannian metric tensor g and Ki/rf 
as an appropriate discrete version of Gaussian curvature K, we [9] defined a combinatorial 
Ricci flow as 

^ = {Rav - Ri)g^, (2.2) 

where gi = rf and Rav = 27rx(M)/||r|p. We further proved that, for surfaces with 
non-positive Euler characteristic, the flow (j2.3p converges iff. there exists a constant 2- 
curvature metric. For a-curvature, we proved similar results for surface M with axiM) < 
0 by introducing a modified flow 

^ = {Ra,av - Ra,i)ri, (2.3) 

at 

where Ra,av = 2vrx(M)/||r||“. However, there are no similar results for surface M with 
axiM) > 0. Even worse, we don’t know whether the solution of flow (j2.2l) exists for all 
time t G (— 00 , -|-oo). For this reason, we introduce another flow here, the solution of which 
always exists for all time t G R. 

Set Ui = In rj. 

Definition 2.2. Given (M, T, 4>) with circle packing metric r, denote Sa = 27rx(M)/||r||“. 
The ath order combinatorial Ricci (Yamabe) flow (“a-flow” for short) is 

^ = Savf - Ki. (2.4) 

Remark 2. Chow and Luo’s normalized Ricci flow (|1.4p is in fact the 0-flow defined above 
with a = 0. Hence a-flow (|2.4p may be considered as a different normalization of Chow 
and Luo’s flow (|1.3p . 

Formally, a-flow seems plainer and simpler than flow (j2.3p . Furthermore, we have 
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Proposition 2.3. The solution to discrete Yamabe flow (12.41) exists for all time t E H. 

Proof. Notice that |sar"| < 2'k\x{M)\, (2 — d)^ < Ki < 27r, where d = max di 

l<i<N 

and di is the degree of vertex i. Hence |s„r" — Ki\ is uniformly bounded by a constant 
depending only on the topology and the triangulation. Then by the extension theorem for 
solution in ODE theory, the solution to (12.41) exists for all t E (—oo, +oo). □ 

Proposition 2.4. a-flow (12.411 is a negative gradient flow. Moreover, a=iD(t) preserves 
to be a constant along a-flow. 

Proof. It is remarkable that Chow and Luo [3] introduced a functional 

ru _ 

{Kj - Kgy) duj, 

which is called discrete Ricci potential, to prove the convergence of their Ricci flow (11.41) . 
With a slight modification, we define the a-order discrete Ricci potential (“a-potential” 
for short) as 

F{u) = / ~ 

JUQ ^ 

where uq E is arbitrary selected. The a-potential is well defined since Y^{Ki — Sarf)dui 
is a closed differential form. Thus a-flow is in fact ii = —V^E. 

Furthermore, — Ki) = 0 implies that Y^ii^) hence 0^(0 are 

invariant along the a-flow. □ 

For this reason, we always assume the initial circle packing metric r(0) belongs to the 
hypersurface r* = 1 in the following. 

Proposition 2.5. If the solution to a-flow ()2.4I1 converges, then there exists at least a 
constant a-curvature metric. 

Proof. Notice that constant a-curvature metrics are exactly the critical points of 
ODE system (|2.4p . If the flow converges, it converges to its critical points. This implies 
the existence of constant a-curvature metric. □ 

Discrete Laplace operators are closely related to discrete curvature flows. Denote 

, 8{Ku---,Kn) 

as the Jacobian matrix of the curvature map. It’s very interesting that L : CiV) —?> Ciy) 
with / 1 -^ Lf can be considered as the prototype of a type of discrete Laplace operator 
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Him 019], which comes from the dual structure of circle patterns. Chow and Luo proved 
that L is positive semi-definite everywhere in lR>o- Moreover, rank{L) = — 1 and the 

kernel of L is t(l, • • • , 1)^, t E IR. Other properties of L is exhibited in Lemma 14.41 The 
authors once defined [9] a two dimensional a-order combinatorial Laplacian as 



Recall that we have considered each / E Ciy) as a column vector, hence the two dimen¬ 
sional a-Laplacian (12.7p can be written in a matrix form, 

( 2 . 8 ) 

with Aq,/ = —for each / E C'(R), where S = diag{^ri, - ■ ■ ,rAr}. It’s interesting 
that we can consider rf as an analogy of dfi, where d/x is the area element in smooth case. 
We can define a a order inner product (•, •)q, on (M, T, <h) with circle packing metric r by 

N 

{f,h)^ = J2fihir? = h^^^f (2.9) 

i=l 

for real functions f,h £ C{V). Then the a-Laplacian Aq : C{V) — C{V) is self-adjoint 
since 

— (/, Aq,/i)q, 

for any f,h£ C{V). 

Denote the first positive eigenvalue of — Aq as Ai(—A„). The following theorem shows 
that the first positive eigenvalue of a order combinatorial Laplace operator is closely 
related to the behavior of discrete a-flow. 

Theorem 2.6. Assuming Ai(— Aq,) > asa = at all r E lR>o- Then the a-flow 

()2.4I) converges iff. there exists at least a constant a-curvature metric on (M, T, 4>). 

Proof. We give a direct and self-contained proof here. We just need to prove the 
“if’ part. Set = S~' 2 LS“' 2 , then 

Aa = —S 2 AqS 2 j 

which implies that Ai(—A q,) = Ai(Aq). Assuming r* is a constant a-curvature metric. 
Scaling r* to any tr* with t > 0, the corresponding a-curvature is still a constant. Hence 
we may suppose r* belongs to the hypersurface 0^ = 1- Denote u* as the rt-coordinate 
of r*, and a-potential as 

^(^)= / ( 2 . 10 ) 

i=l 
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Denote r“ = (rf, • • • , then it’s easy to calculate 


HesSuF = L — asa 




oe 


= S2 





a 

OL 


OL 


S2. 


( 2 . 11 ) 


Choose an orthonormal matrix P such that P^AaP = diag{0, Xi{Aa), ■ ■ ■ , A 7 v-i(Aq,)}. 
Suppose P = (eo, ei, • • • , cty-i), where e* is the {i + l)-column of P. Then A^eo = 0 and 


AaCi = AjCj, 1 < i < — 1, which implies cq = 

Hence (l - )eo = 0 and (/ - )& 


= e,: 


1^11 1 — 

\r 2 II and r 2 T Cj, 1 < i < — 1. 

1 < * < A^ — Ij which implies 


HesSuF = S 2 Pdiag{0, Ai(Aq) — aSa, • • • , XN-i{Aa) — aSajS 2 


If Ai(Aq,) > aSa = a^Tpri^, then HessF > 0, rank{F) = N — 1. Using Lemma 
14.11 Lemma 14.21 and Lemma 14.31 we know that F is proper on ^, where '2^ = {n G 
= 0}. Moreover, lim F{u) = +00 and u* is the unique zero point of 

, u^oo 

VP which is also the unique minimum point of F on Let ip{t) = F{u{t)), then 
= — Yli{Ki — Savf)^ = — ||VP|p < 0, which implies that u{t) lies in a compact 
region of ^. We can further get = 2{K — So,r°‘)^HessF{K — So,r^) > 0. The fact 

that if' < 0, (/9" > 0 and (p is bounded below implies (/9'(+oo) = 0. Combining u* is the 
unique zero point of VP and {u{t)} CC ^, we know u{t) — u* exponentially fast. □ 


If ax(M) < 0, we always have Ai(—A q,) > 0 > aSa = thus we have 

Theorem 2.7. Suppose (M, T, 4>) is a weighted triangulated surface with ax(M) < 0. 
Then the solution to the a-flow (Hai) converges if and only if there exists a constant a- 
curvature metric r*. Furthermore, if the solution converges, it converges exponentially 
fast to the metric of constant curvature. 

Remark 3. When a = 0, above theorem is obtained in [3]. 

Next we want to study the combinatorial and topological conditions for the existence 
of constant a-curvature metric. For any proper subset I C V, denote 

‘3^1 = {x ^ (tt - $(e)) + 27rx(P/)} 

iGl (e^v)£Lk{I) 

and JHgb — {x G R^| = 27rx(^)}- The authors once proved 

Proposition 2.8. ([9], Theorem 2.34) Given a weighted triangulated surface (M, T, $). 
Consider Thurston’s circle packing metric and the a-curvature. When a > 0 and x(Af) < 
0, the existence of constant a-curvature metric is equivalent to ^)CR<o 
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(f); When a < 0 and xi^) > 0, the existence of constant a-curvature metric is equivalent 
to J(^GB n n lil >0 ^ 4>'i When x{M) = 0) the existence of constant a-curvature 

metric is equivalent to H '3^i) n {0} / 4>. 

Similar to Thurston’s condition o, conditions in Proposition (12.Sp also show that 
the combinatorial structure of the triangulation and the topology of surface, which have no 
relation with circle packing metrics, contains a-curvature information surprisingly. Using 
ProDosition l2.8l we can derive a combinatorial-topological-metric condition which contains 
Thurston’s condition (11.21) as a special case. 

Theorem 2.9. Suppose (M, T, <h) is a weighted triangulated surface with ax{M) < 0. 
Then there exists a constant a-curvature metric iff. there exists a circle packing metric r* 
such that for any nonempty proper subset I of vertices V, 

27rx(M) ^^y^j >- ^ (,r-cJ>(e)) + 27rx(T,). (2.12) 

{e,v)eLk{I) 

Proof. By Lemma 14.51 for any circle packing metric r and any proper subset I dV, 
the classical combinatorial Gauss curvature K satisfies 

- ^(e)) + 2vrx(T/). (2.13) 

iGl (e^v)GLk{I) 

If there exists a constant a-curvature metric r*, substituting K* = K{r*) = {K*, • • • , K'^) 
into (|2.13l) . where K* = we get (|2.12l) . On the other hand, it’s easy to see, 

• when a > 0 and x{^) < 0; (|2.12l) implies that J^gb O (n^^jcy ^) n ^ 0; 

• when a < 0 and x{M) > 0, (|2.12l) implies that J^gb O (n^^jcy ^) n R^q ^ cj)] 

• when x{M) = 0) (I2.12h implies that 0 G J^gb O (O^^/cy 

Proposition 12.81 shows that conditions in above three cases all implies the existence of 
constant a-curvature metrics. For a = 0 case, (12.121) is in fact Thurston’s condition (II.2|) . 
Thus we finish the proof. □ 


Remark 4. The proof of Proposition 12.81 and hence Theorem 12.91 deeply rely on deriving 
a discrete maximum principle for the flow (12.31) . We want to know if there are more direct 
proofs without using discrete maximum principle for the flow (12.3|) . 

As to the uniqueness of constant a-curvature metric, we restate Theorem 2.33 [9] here 
for completeness. 

Theorem 2.10. Suppose (M, T, $) is a weighted triangulated surface with ax(M) < 0, 
then the constant a-curvature metric is unique if it exists. Specificly, if ax(M) = 0, then 
there exists at most one constant a-curvature metric up to scaling. If ax(Af) < 0, then 
for any c*, there exists at most one metric with a-curvature Ra,i = c*. 
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When ax{M) > 0, such as a = 2 and M is a sphere, Example 2 in [9] shows that 
the conclusions in Theorem 12.71 are not true. For the tetrahedron triangulation of the 
sphere, if the initial metric r(0) is close enough to r* = (1, • • • , 1)^, then the solution to 
flow (j2.4p converges to r* when t —)• —oo. However, the limit behavior of r(t) depends on 
the selection of initial metric r(0). Indeed, there exists r(0) such that the solution r{t) 
diverges to oo either t tends to +oo or —oo. Example 3 in [9] shows that the constant 
a-curvature metric is not unique. For the existence of constant a-curvature, we have 

Corollary 2.11. Suppose (M, T, <I>) is a weighted triangulated surface with ax(M) > 0. 

(1) There exists a constant a-curvature metric. 

(2) There exists a circle packing metric r* such that for any nonempty proper subset I of 
vertices V, 

>- {^-He)) + 27rx{Fj). (2.14) 

{e,v)eLk(I) 

(3) When a > 0 and x(M) > 0, then JFgb H (n^^/cy '3^i) H / (p; When a < 0 and 
x(M) < 0, then jecB n (n^^/cy n E^g ^ 4>. 

Then (1) implies (2) which implies (3). 

Proof. For (1) ^ (2), using formula ()2.13p . For (2) ^ (3), it’s obviously. □ 


2vrx(M) 


E 


i&I 


Definition 2.12. Suppose (M, T, ‘h) is a weighted triangulated surface with circle packing 
metric r, R G C{V) is a function defined on M. The modified a-flow with respect to R is 
defined to be 

^ = R,rf-K,. (2.15) 

R is called admissible if there is a circle packing metric r with curvature R. The 
modified a-flow can be used to study prescribing curvature problem. On one hand, if the 
solution to the modified a-flow (I2.15P converges, then R is admissible. On the other hand, 
we have 

Theorem 2.13. Suppose {M,T, 4>) is a weighted triangulated surface and R E C{V) is a 
function defined on M. If aRi < 0 for all i, but not identically zero, and R is admissible by 
a metric r. Then r is the unique metric in E>g such that it’s a-curvature is R. Moreover, 
the solution to the modified flow (|2.15p converges exponentially fast to r. 

Proof. The first part is obviously, and R is admissible by metric r(-|-oo). For the 
second part, for given function R E C'(F), we can introduce the following modified a- 
potential 

ru ^ 

(2.16) 

Juo j=l 
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It is easy to check that the modified a-potential F is well-defined. Furthermore, by direct 
calculation, we have 


HesSuF = L — T,2 


aRi 


E2. 


oRn / 


It is easy to check that, if aRi < 0 for i = 1, • • • , and not identically zero, HesSuF is 
positive definite. By Lemma [l]3l = {Ki — i?irf, • • • , iLjv — R^rff)'^ is an injective 

map from u G to R^. Hence r is the unique zero point of V„F. This fact implies 
that r is the unique metric in R^q such that it’s curvature is R. By Lemma 14.21 we know 
that F is proper and lim F{u) = -|-oo. Furthermore, ^F{u{t)) = — ~ < 0 

implies that the solution of (|2.15l) lies in a compact region. The following of the proof is 
the same as that of Theorem 12.61 so we omit it here. □ 


Remark 5. The second part of Theorem 12.131 implies that ax(M) < 0. If aRi = 0 for 
all i, then the corresponding prescribing curvature problem is already solved in Theorem 
O In this case, the metric r is not unique. However, it’s unique up to scaling. This is 
slightly different from Theorem 12.131 


In the following of this section, we consider more general “area element” Ai. It’s 
interesting to define the “H-curvature” as Ri = KijAi, where > 0 is a function of circle 
packing metric r G R>o- We can consider the following discrete flow 


Kii) 


27rx(M) 


■^i Fj, 


(2.17) 


which is called “H-flow” for short. This generalized H-flow can be used to evolve a metric 
to a metric with constant H-curvature, i.e. a metric r satisfying Ki = sAi for all i G H, 
where s = It’s easy to see the solution to H-flow always exists for all time 

t G (—oo, -|-oo). So this flow seems simpler than other flows such as iii = s — Ri, which is 
an H-generalization of flow (12.3p . Furthermore, if the solution r{t) to H-flow converges to 
r(+oo) G R> 0 ) then r(+oo) has constant H-curvature. It’s very interesting that we can 
select Ai as the real area instead of the area of disk packed at i, hence 


N 

Ai = Area{M, T, <h, r) (2-18) 

i=l 


is necessary, where Area{M, T, r) is the total real area of a weighted triangulated surface 
(M, T, ‘h) with a fixed circle packing metric r. Then it’s easy to see the following three 
selections of Ai all satisfy (j2.18ll . 
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Example 1. Select Ai = Area(Aijk)/3, where the sum is taken over all the 

/\ijk£F 

triangles with i as one of its vertices. 

Example 2. Consider the dual structure determined by Thurston’s circle patterns. For 
any face {ijk} G F, denote Ci,Cj,Ck as the closed disks centered at i, j and k so that 
their radii are r,, rj and r^. They both intersect with each other. Let Ci^ Cj, Ck be the 
geodesic lines passing through the pairs of the intersection points of {(7^,(7^}, {(7^,(7*}, 
{Ci,Cj}. These three lines £j, Ck must intersect in a common point Oijk- Connect 
Oijk and Oiji whenever triangles {ijk} and {ijl} share a common edge {ij} G E. Thus 
we get a dual graph. For more details see mill [12]. Select Ai = Area{Di), where Di is 
the dual 2-cell of i. 


Example 3. Inspired by |16] . we can select Ai = Area{Vi), where Vi is the Voronoi dual 
2-cell of i in the Delaunay triangulation. 


Under the assumption (|2.18p . the constant A-curvature is ^ r) • 

approximates a smooth Riemannian surface {M,g), then Area{M,F,^,r) approximates 
Area{M,g). Hence the constant H-curvature # r) approximates the smooth aver¬ 


age curvature 


This fact inspires us to consider the following problem. 


Conjecture 1. Fix a smooth Riemannian surface {M,g). Suppose {Mn,Tn,^n,'f’n) is 
a sequence of weighted triangulation of M with initial circle packing metric r„. is 
different with {M,g) as metric space, although they are topologically equal. For each n, 
one can evolve H-flow (j2.17p and get a solution t G [0,-|-oo). Meanwhile, one can 

evolve {M,g) by smooth Ricci flow and derive a solution g(t), t G [0, -|-oo). Assuming the 
initial Gromov-Hausdorff distance between (M, g) and (M„, Tn, r^) tends to zero, then 
^nit) —^ g(t) as n ^ -boo. 


3 a order combinatorial flow in three dimension 


Suppose M is a 3-dimensional compact manifold with a triangulation T = {V,E,F,T}, 
where the symbols V, E, F, T represent the sets of vertices, edges, faces and tetrahedrons 
respectively. A sphere packing metric is a map r : U ^ (0, -boo) such that the length 
between vertices i and j is kj = r* -b rj for each edge {i, j} G E, and the lengths 
hj:hk:kh ^jk,lji,hi determines a Euclidean tetrahedron for each tetrahedron {i, j, k, 1} G T. 
Glickenstein pointed out [inj that a tetrahedron {i,j,k,l} G T generated by four positive 
radii rj, rj ,rk,ri can be realized as a Euclidean tetrahedron if and only if 


Qijkl 



> 0 . 


(3.1) 
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Thus the space of admissible Euclidean sphere packing metrics is 


TIt = { r G I Qijki > 0, G T } . 

Cooper and Rivin [4] called the tetrahedrons generated in this way conformal and proved 
that a tetrahedron is conformal if and only if there exists a unique sphere tangent to all 
of the edges of the tetrahedron. Moreover, the point of tangency with the edge {*,j} is 
of distance r* to Uj. They further proved that OJIt- is a simply connected open subset of 
IR^g, but not convex. 

For a triangulated 3-manifold (M, T) with sphere packing metric r, there is also the 
notion of combinatorial scalar curvature. Cooper and Rivin [1] defined combinatorial 
scalar curvature Ki at a vertex i as angle deficit of solid angles 


dvr ^ ^ Olijkh 

{i,j,k,i}eT 


(3.2) 


where aijki is the solid angle at the vertex i of the Euclidean tetrahedron {i,j,k,l} G T 
and the sum is taken over all tetrahedrons with i as one of its vertices. For this curvature, 
Glickenstein [TO] first defined a combinatorial Yamabe flow 

§ = (3.3) 


and give some very interesting and inspiring results. 

Similar to the two dimensional case, Cooper and Rivin’s definition of combinatorial 
scalar curvature Ki is scaling invariant, which is not so satisfactory. The authors [9] once 
dehned a new combinatorial scalar curvatures as Ri = Ki/rf on 3-dimensional triangulated 
manifold (M, T) with sphere packing metric r. Consider rf as the analogue of the smooth 
Riemannian metric. If = erf for some positive constant c, we have Ri = c~^Ri. This 
is similar to the transformation of scalar curvature in smooth case under scaling. For 
this type of combinatorial scalar curvature, the authors defined a combinatorial Yamabe 
functional 


Q{r) = 

Vs 


Eh 


r G fUIri 


(3.4) 


(E£i 

and proposed to study the corresponding constant curvature problem which is called com¬ 
binatorial Yamabe problem. For this, the authors defined a new discrete Yamabe flow 


with normalization 


dgi 

dt 


— Ridi ) 


dgi 

dt 


{Rav Ri)gij 


(3.5) 
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(3.6) 



where gi = rf and Rav = —3 is the average of the combinatorial scalar curvature. 

Constant i?-curvature metric means that Ri =constant, which implies K = Rayr^- The 
authors [7] once dehned the so called discrete quasi-Einstein metric satisfying K = \r, 
which is similar to constant i?-curvature metric. Motivated by these phenomena, we can 
generalize these properties to a order combinatorial scalar curvature (“a-curvature” for 
short). 


Definition 3.1. For a triangulated 3-manifold (M, T) with sphere packing metric r, the 
a-curvature at the vertex i is dehned as 




(3.7) 


for any a € R, where Ki is given by (13.20 . 


The study of smooth Einstein-Hilbert functional has a long history. For piecewise 
linear metric case, Regge m hrst give a discretization of this functional. For sphere 
packing metric case, the Einstein-Hilbert-Regge functional S = YliLi introduced 

by Cooper and Rivin in [3]. 


Definition 3.2. Suppose {M,'T) is a triangulated 3-manifold with a hxed triangulation 
T. For any a G R, —1, the a order combinatorial Yamabe functional (“a-functional” 
for short) is dehned as 




Qa{r) = 


rllo+i 


:e£i 


1 

+1 


r G SOIt". 


(3.8) 


The a order combinatorial Yamabe invariant with respect to T is dehned as 


Ymj = inf Qair), 

r&Kq- 

while the a order combinatorial Yamabe constant of M is dehned as Ym = sup inf Qair). 

r r€Wr 

When a > 0, then |Qa(r)| < ||R||^,i. When —1 < a < 0, then |Qo(r)| < Ei 177*1. 

a 

Hence the a order combinatorial Yamabe invariant Ym,t is well dehned when a > — 1. 
For a > 0 case, Ym,t attains the minimum value — ||i7||i, 1 at a metric r* if and only if 

’ OL 

r* is a constant a-curvature metric with s* < 0 . As noted in [9], the admissible sphere 
packing metric space TUt- for a given triangulated manifold (M, T) may be considered as 
the combinatorial conformal class for (M, T), which is an analogue of the conformal class 
[g'o] of a Riemannian manifold (M, g'o)- Denote Sq, = 5/||r||“^}. Then we have 

^rQa = II II . (3.9) 
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Hence r is a constant a-curvature metric iff. it is a critical point of a-functional Qa{r). We 
raise the following discrete a-Yamabe problem on 3-dimensional triangulated manifold. 

The Combinatorial a-Yamabe Problem. Given a 3-dimensional manifold M with 
triangulation T, find a sphere packing metric with constant combinatorial a-curvature in 
the combinatorial conformal class SIRt-. 

It’s easy to see, if r is a constant a-metric with K = Ar", then A = Sq, = 5/||r||"^j;. To 
study the combinatorial a-Yamabe problem, we introduce 3-dimensional a-Yamabe flow. 

Definition 3.3. Given a triangulated 3-manifold (M, T) with sphere packing metric r. 
For any a G IR, the a order combinatorial Yamabe flow (“a-flow” for short) is 

^ = SaVf - Ki. (3.10) 

at 

Remark 6. The flow r = \r — K (A = si) introduced by the authors in [7] is in fact the 
a-flow defined above with a = 1. Hence a-flow (j3.10p may be considered as a different 
normalization of the flow r = Xr — K. 

Along the 3-dimensional a-flow, ||r(t)||| = is invariant. Hence we always 

assume r(0) G in the following. It’s easy to see that, if the solution of (I3.10p converges 
to a metric r(oo), then r(oo) is a metric with constant a-curvature. It’s interesting that 
almost all 3-dimensional results in [3 [9] are still true for a-flow (I3.10p . We just state some 
of them here. 

Theorem 3.4. If the solution of (|3.1Up lies in a compact region in 9Jl7-nS'^“^, then there 
exists at least one sphere packing metric with constant a-curvature on (M, T). □ 


Definition 3.5. Given a triangulated 3-manifold (M, T). For any a G R, the a order 
combinatorial Laplacian (“a-Laplacian” for short) : C(V) —)• C(V) is defined as 


Aafi 


1 


E( 


dKj 

drj 


;■)(/. 


fi) 


(3.11) 


for / G C{V). 


This definition of a-Laplacian is a generalization of a = 2 case, which was carefully 
studied by the authors in [^. Similar to the two dimensional a-Laplacian, the three 
dimensional a-Laplacian (j3.11l) can also be written in a matrix form. 




(3.12) 
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(3.13) 


with Aaf = —S “AS/ for each / G C{V), where S = diag{ri, ■ ■ ■ ,rAr} and 


A = HesSrS 


,Kn) 

d{ri, ■■■ ,rN) 


It was proved diiHiiinKii] that A is positive semi-dehnite with rank N —1 and the kernel 
of A is the linear space spanned by the vector r (see Lemma 14.61) . 


Set rj(r) = Savf — Ki, 1 < i < N. Then the a-flow (13.10^ can be written as r = r(r), 
which is an autonomy ODE system. Differentiate T, we get 


Dr(r) 


-A T CXSq: 


^diag{r^ \--- 


,r 


OL—l 

N 


} 


|| 7-||“+1 J 
II' lla+l / 


{K- 




IQ+l 

la+l 


If r* G TUt" is a sphere packing metric with constant a-curvature, then 


af a\T \ \ 

diag{r°-\--- ^ ^ . (3.14) 

fL +1 

Proposition 3.6. Given suppose r* is a constant a-curvature metric. If the 

first positive eigenvalue of — Aq, at r* satisfies 



Ai(-A„) > as* 


a5* 


II *||“+i 

II' llo+l 


(3.15) 


then —DT\r* is positive semi-definite with rank N — 1 and kernel {tr*\t G IR,}. 
Proof. Denote A = AS^“. Then 


_ . 1 + q: ~ l + a 

-A„ = S-“AS = S-—AS— 


which implies that 

Ai(-A„) = Ai(A). 


Choose a matrix P G 0{N) such that P^AP = diag{0, Xi{A), ■ ■ ■ ,AAr_i(A)}. Suppose 
P = (eo,ei,--- ,eAr_i), where e* is the {i + l)-column of P. Then Acq = 0 and Ae* = 


Q + 1 


AjCj, 1 < f < A — 1, which implies eg = r 2 


I 

r 2 


I- 


, S+i.T’ 

—iirir+ 1 — 

Ip lla+l 


eg = 0 and I — 


Q +1 ^ 0+1 

-Ilrll“+1- 

Ip lla+l 


a + l 

and r^~ T e^, 1 < i < A — 1. Hence 
= ej, l<z<A—1. Furthermore, 


—Dr|r* = S 2 Pdia^f |o, Ai(A) — as* ,•••, AAr_i(A) — as* I S 2 p^. 
Hence the conclusion is derived. 


□ 
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Theorem 3.7. The constant a-curvature metrics satisfying Ai(—A q,) > asa are isolated 
in Specifically, the constant a-curvature metrics with aSa < 0 are isolated. 

Proof. Consider the map F : WIt —>■ IR'^, r i—)• r(r). It is easy to see that the zero 
point of r corresponds to the metric with constant a-curvature. By (I3.14p and Proposi¬ 
tion [21S1 if the conditions in this theorem are satisfied, then DT\r*, the Jacobian of T at 
the constant a-curvature metric r*, is negative semi-definite with rank N — 1 and kernel 
{tr*\t E R}. Notice that the kernel is the normal of Restricted to DT is 

negative definite and then nondegenerate, which implies the conclusions. □ 

Proposition 13.61 shows that constant a-curvature metric r* with Ai(—A*) > as* is a 
asymptotically stable point of the a-flow. Hence we have 

Theorem 3.8. Given a triangulated manifold (M^,T). Suppose r* E is a constant 

a-curvature metric satisfying Ai(—A*) > as*, or more specifically, r* E is a con¬ 

stant a-curvature metric with as* < 0. If ||r(0) — r*|| is small enough, then the solution 
of a-flow (I3.10p exists for all time and converges to r*. □ 


Let’s look at what happens when a = 0. The 0-curvature is just Cooper and Rivin’s 
scalar curvature. By Theorem 13.71 and Theorem 13.81 constant curvature metrics r* are 
always isolated in and are always local attractors of the ODE system (I3.10p . 

Remark 7. a-flow ()3.10p is not a negative gradient flow, since d{sarf)/drj ^ d{sar‘j)/dri. 
We could get similar results by considering the negative gradient flow of a-functional (j3.8p , 
i.e. ^ 

r = -VrQa = - {Sarf - K). (3.16) 

For this flow, Theorem 13.41 and Theorem 13.81 are still true. 

4 Some useful lemmas 

For reader’s convenience, we list some lemmas in this section which are used in the proof 
of our main results. Some of them are proved in detail while some of them we just give 
related references. 

Lemma 4.1. Given a function F E C'^(R'^). Assuming Hess{F) > 0, rank{Hess{F)) = 
A — 1, VF has at least one zero point and F{u + t(l, • • • , 1)"^) = F{u) for any u E 
and t E R. Denote '^ = {u E R^| Yli — 0}- Then the Hessian of F\^ (considered as a 
function of A — 1 variables) is positive definite. 
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Proof. Chow and Luo [3] stated similar results for a-potential F with 0 = 0. Since 
this lemma is very important, we give a direct and rigorous proof here for completeness. 
Set 7 = (1, • • • , l)^/\/]V. Choose A E 0{N), such that ^47 = (0, • • • ,0,1)^, meanwhile, 
A transforms ^ to {C E R^\Cn = 0}. Define ,Cn-i) = f{A'^{Ci,--- , Cat-i, 0)'^), 

we can finish the proof by showing that Hess{g) is positive definite. Partition A into two 
blocks, A = 7 ]. Write L = Hess{F) for short, then Hess{g) = S^LS. Notice 

that, F{u + t'j) = F{u) implies L 7 = 0. Therefore we have 


A^LA 



L[S,j] 


'S^LS 


'S'^LS o' 

7^L7 


0 0 


which implies that S'^LS is positive semi-definite. Furthermore, rank{S'^LS) = N — 1. 
Hence Hess{g) = S'^LS is positive definite, which implies that F\^ is strictly convex. □ 


Lemma 4.2. Given a function F E C'^(1R'^) with Hess{F) > 0. Assuming VF has at 
least one zero point, then lim F(x) = -|-oo. Moreover, F is proper. 

||a:||—)-oo 

Proof. Without loss of generality, assume VF(0) = 0. First, we consider the case 
N = 1. So we have F : lH —)• IR with F" > 0 and F^(0) = 0. Since F" > 0 and F^(0) = 0, 
then F'(l) = a > 0 and F'(—1) = —6 < 0. Again, by F" > 0, F'{x) > F'{1) = a > 0 for 
X > 1, and F'{x) < F'{—1) = —b. Integrating this gives 

F{x) > a(x — 1) + F(l) for X > 1 
F{x) > b{—l — x) -|- F(—1) for x < —1. 

That means F{x) > min{a, b}\x\ +C, where C = min{—a + F(l), — 6 -|-F(—1)}. Obviously, 
lim F{x) = +00 and F is proper. 

|a;|—>-oo 

For > 1, for any oj E we consider the ray {tw}o<t<cxD C Let = 

F(tco), then /(^(O) = 0 and f"{t) = HessF{u},u}) > 0, then as in the dimension N = 1 
case, there exists constant a^j = /^(l) = VF(a;) ■ ui >0 and such that 

F{tu}) = fuj{t) > ttujt — acj + F{lo), for t > 1. 

By the arguments above, we have proved a;-VF(w) > 0 for any a; E by compactness 

of and F is C^, 

A := inf oj ■ \/F(u}) > 0. 

Let B := min F(uj), then we have 

F{x) > A||x|| — A + B, for ||x|| > 1. 
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It implies lim F{x) = +oo and F is proper. □ 

||x||—>-oo 

Lemma 4.3. ([3]) Suppose H C is convex, the function /i : —>■ IR is strictly convex, 

then the map \7h : Q ^ IR'^ is injective. 

Lemma 4.4. ([3], Proposition 3.9.) Suppose (M, T, is a weighted triangulated surface. 
r = (ri, • • • , is circle packing metric, while K = {Ki, • • • , is classical discrete 

Gauss curvature. Then L = is positive semi-definite with rank N-1 and kernel 

{t(l, • • • , l)^|t G H}. Moreover, > 0, < 0 for i ~ j and = 0 for others. 

Lemma 4.5. ([HI [l5l |3]) Given a weighted triangulated surface (M, T, <h). Gonsider 
Thurston’s circle packing metric r and classical discrete Gauss curvature K. For any proper 
subset I C V, denote ^ = {x G > -T.(e,v)GLk{i)i'^ “ ^(e)) + 27rx(F/)}, 

^GB — {x G IR'^I = 27rx(-Y)}. Then the space of all admissible iF-curvature is 

{K = Kir)\r G IR^ol = <^GB G ^i). 

Lemma 4.6. (jH [HI [HI [H] ) Suppose {M,T) is a triangulated 3-manifold with sphere 
packing metric r, 5 = H is the Einstein-Hilbert-Regge functional. Then we have 

VrS = K. (4.1) 

dKi \ 

9rjv 


. . dKi, 

drj^ / 

kernel of A is the linear space 
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If we set 


/ 

dr\_ 


A = HesSrS = 


,Kn) 

d{ri, ■■■ ,rN) 


V 


OKn 

dr I 


then A is positive semi-definite with rank N — 1 and the 
spanned by the vector r. 
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